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I. INTRODUCTION 



The study of parity violating (PV) and time reversal invariance violating (TRIV) effects in 
low energy physics is very important for understanding main features of the Standard model 
and for a search for new physics. During the past 50 years, many calculations of different 
PV and TRIV effects in nuclear physics have been done. The methods of calculations 
of PV and TRIV effects are very similar, therefore, the reliability of prediction of TRIV 



_PV effects. Unfortunately, in the 
and references therein) that the 



effects can be justified by the accuracy of calculation o: 
last few years, it has become clear (see, for example 
traditional Desplanques, Donoghue, and Holstein (DDH) jsj method for the calculation of 
PV effects cannot reliably describe the available experimental data. It could be blamed 
on "wrong" experimental data; however, it may be that DDH approach is not adequate 
for the description of the set of precise experimental data because it is based on a number 
of models and assumptions. Recently, a new approach based on the effective field theory 



(EFT) 



las been introduced as a model independent parameterization of PV effects (see, 
,y] and references therein), and some calculations for two-body systems have been 



papers 

done Q|. The power of the EFT approach could be utilized if we can analyze a large enough 
number of PV effects to be able to constrain all free parameters of the theory, which are 
usually called low energy constants (LEG), to guarantee the adequate description of the 
strong interaction hadronic part of symmetry violating observables. Then, if discrepancies 
between experimental data and EFT calculations will persist, it will be a clear indication 
that the problems are related to weak interactions in nuclei and probably to a manifestation 
of new physics. 

Unfortunately, the number of experimentally measured (and independent in terms of 
unknown constants) PV effects in two body systems is not enough to practically constrain 
all LECs. Therefore, one has to include into analysis few-body systems and even heavier 
nuclei, which are actually preferable from an experimental point of view, because usually, the 
measured effects in nuclei are much larger than in nucleon-nucleon scattering due to nuclear 



enhancement factors 



To verify the applicability of the EFT approach for calculations 
of symmetry violating effects in nuclear reactions, it is natural to start from a scattering 
problem in three-body systems, and to develop a regular and self-consistent approach for 
calculation of symmetry violating amplitudes in a few-body systems, which later could be 



extended to many-body systems. From this point of view, we present an analysis of different 
methods of calculation of symmetry violating effects in few-body nuclear systems. Taking 
into account the similarity of techniques for calculation PV and TRIV effects, we will refer 
further only to PV effects assuming that the same could be applied for calculation of TRIV 
effects. 



II. METHODS OF CALCULATIONS 



Because of small value of weak coupling constants, PV effects in nuclei can be calcu- 
lated in the first order of perturbation theory and represented as a sum of terms multiplied 
by corresponding weak interaction related constants. For example, in the DDH approach, 
nuclear PV effect is represented as a linear superposition of terms multiplied by correspond- 
ing weak nucleon-meson coupling constants. In EFT approach, one can expect to obtain 
a description of nuclear PV effects as a linear superposition of terms multiplied by small 
numbers of LECs. To obtain PV amplitudes, we may introduce the weak potential derived 
from EFT Lagrangian or directly sum Feynman diagrams from EFT Lagrangian. Although 
the derivation of PV EFT potential would not be unique, the calculation of two body PV 
amplitudes can be well defined. Thus, let us call the first approach the "hybrid" method and 
the other one the "true" EFT method. For all these cases with potentials, the important 
issue is how to calculate nuclear wave functions: to use Schrodinger equations or Faddeev- 
type few body equations. It should be noted that, due to small value of weak interactions. 
Distorted Wave Born Approximation (DWBA) is a standard method for the calculation of 
PV amplitudes using wave-function obtained both from Schrodinger and Faddeev equations. 
Therefore, one can classify all methods of calculations first by the method of description of 
two-body weak interactions: DDH or EFT methods. The DDH potential can be used to 
calculate amplitudes with wave functions obtained from Schrodinger or Faddeev equations, 
which leads to "nuclear reaction" or "few-body equations" approaches. 



A. "True" EFT approach 

To relate PV reaction matrices (amplitudes) for neutron nuclei interactions with two- 
nucleon weak reaction matrices, one can use AGS-type jlo| of Faddeev equations ^Jj] written 



in terms of transition operators Upa'- 

= ^al3GQ^ + ^ S^ij^GqU^^, (l) 
7 

where upper idex s indicates strong interactions, Sajs = l — Sa/s, is a two-particle transition 
operator in three-particle space, a,/3, . . . are channel index, and 

Goiz) = -i- (2) 
z — rin 



is the resolvent operator for free motion. Then, as it was shown in paper [12|, the PV 
transition operator can be expressed in terms of PV two-particle scattering operator 
using the integral equation 

^Ja = '^7/3^7^of^7« + ^ ^-fpf^GoU^^, (3) 

7 7 

where two-body scattering operators are written as a sum of PC and PV (indicated by w 
for weak interactions) parts: 

t^ = t; + t!;. (4) 

It should be noted that the first term of the integral equation([3]) has exactly the same kernel 
as the equation ([T]) which describes the process with strong interactions only. The second 
term of Eq.(l3]) does not depend on f/J^ and, therefore, it is a free term in the integral 
equations. One can see that the first (integral) term includes strong two-body transition 
operators but the second one (free term) contains a direct contribution to UJ^^ from weak 
two-body transition operators. Therefore, Eq.(l3]) gives us a framework for a calculation of 
symmetry violating amplitudes using Faddeev type three-body equations in terms of two- 



body amplitudes, which can be obtained from EFT 
has not been applied to few-body systems yet. 



Unfortunately, this method 



B. Faddeev equations with DDH potentials 



The natural way to calculate PV effects in the DDH formalism for few-body systems 
is to calculate PV amplitudes using DWBA with wave-function obtained from few-body 



equations (see, for example [l3|, ll^). In this case, PV effects can be represented as a 



linear superposition of contributions from different parts of the DDH potential weighted by 



corresponding weak meson-nucleon coupling constants. To see how PV effects are sensitive 
to these DDH couphng constants, one can consider results of calculations of PV neutron spin 
rotation and neutron spin asymmetry P in the propagation of polarized neutrons through 
unpolarized deuteron target. We define the angle of neutron spin rotation per unit length 
of the target sample in terms of elastic scattering amplitudes at zero angle for opposite 
helicities /+ and /_ as 

d<j) 2nN .... 

-r = Re /+-/_ , (5) 

dz p 

where is a number of target nuclei per unit volume and p is a relative neutron momentum. 
The neutron spin asymmetry is equal to the relative difference of total cross sections P = 
— (T_)/ ((T_|_ + for neutrons with opposite hehcities. 
Then, one can calculate these PV effects for three different choices of DDH constants in 
Table 



D which correspond the the DDH "best" values, 4-parameter fit [15|, and 3-parameter 
fit [15|]. As a result of the calculations, PV neutron spin rotation angles have about the 
same value 7.68, —6.82 and —8.91 in units of 10^^ rad-cm^^ for these three choices of weak 
coupling constants, correspondingly. The values of neutron spin asymmetry are 8.99, —2.09, 
and —2.36 in units of 10~^. This shows that different PV parameters have, in general, 
different sensitivity to the DDH coupling constants. Therefore, making a large enough 
number of high accuracy measurements of different PV effects, one can, in theory, constrain 
values of weak coupling constants. However, it is a very difficult task because values of PV 
effects in few-body nuclei usually are very small. 



C. "Hybrid" EFT approach 

The "hybrid" approach uses the same calculation scheme, as it was described in the 
previous section, but instead of using DDH potentials, it uses the potentials derived from 
a particular choice of EFT Lagrangian. However, as it was shown in [l3f|, the potentials 
derived from pionless EFT lagrangian [l| contain contributions from the same operators 
as the DDH potentials [16]. This is also correct for potentials derived from pionful EFT 
Lagrangian [li]. Therefore, all these potentials (DDH, EFT pionless and EFT pionful) can 
be expanded in terms of a set of O-"'' operators as 

vfj = CnOlf, a = DDH or pionless EFT or pionful EFT (6) 

n 

5 



where operators 0\j are defined as products of isospin, spin, and vector operators X -^'^ 

x^l^z[p,^,Un,)U (7) 

witli = (pj— Pj)/2, and c" being expansion parameters. Tlie parameters of this expansion 
or DDH, EFT pionless, and EFT pionful potentials are presented in Table HTl (see paper 
1^ for more details and discussions). One can see that all these potentials (DDH and 
EFT ones) have phenomenological radial dependencies (/^(r), L\{r), Lf^{r), and Hf^{r) ) 
which are related to masses of exchange mesons for the DDH case and to cutoff parameters 
in the case of EFT potentials (for detailed notations and discussions of the properties of 
these radial functions, see |l^). Therefore, the results of calculations of PV effects in the 
"hybrid" approach can be represented as a linear combination of contributions from different 
operators weighted by (unknown) LECs and folded (integrated out) parameters which are 
depended on these radial functions. This is very similar (see Table |Tl]) to the representation 
of PV effects of the DDH- type of calculation. The only practical difference of the "hybrid" 
approach, as compare to the DDH one, is related to the different numbers of operators 
corresponding to different models of EFT. For example, for DDH potentials, it could be up 
to 13 operators, while for considered EFT potentials, it could be 5 or 9 operators for pionless 
and pionful EFTs, correspondingly. 

Despite the similarity of structure between DDH and "hybrid" EFT potentials, DDH 
potential has additional assumptions and constraints on the forms and relations between 
operators. Because the systematic fitting between "hybrid" EFT calculation and experi- 
mental data has not been done yet, we cannot make any prediction using "hybrid" EFT. 
Also the adoption of "hybrid" EFT may bring some undesirable artificial effects. Thus, 
"true" EFT calculation may be an ideal approach to the few-body PV observables. 

D. Nuclear reaction approach 

Since PV effects in a few-body systems are usually too small to expect many experimental 
results and precise calculations of these effects are rather difficult, it is desirable to have a 
method for a reliable estimate of possible observable parameters. It is possible to calculate 
some PV effects in the first order of perturbation theory with wave functions obtained from 

6 



a solution of Schrodinger equation for reliable nuclear models (see, for example Q, Q). 
We consider here a nuclear reaction approach which gives PV amplitudes in DWBA with 
wave functions obtained from a reliable nuclear models with parameters fixed from available 
experimental data. This gives the opportunity to explore many nuclei and observables and to 
find nuclei with reasonably large PV effects for further experiments and detailed calculations. 
To illustrate this approach, let us estimate PV effects in n +^ He — J-^ H + p reaction with 
polarized neutrons. The ^He and ^He systems were subjects of intensive investigations 
for a long time, and as a result, many parameters related to reactions with neutrons and 
protons, as well as to excitation energy levels of these nuclei, have been measured and 
evaluated by a number of different groups. This rather comprehensive data provides the 
opportunity to estimate values of possible PV effects and their dependence on neutron 
ener gy i n n +^ He — t-^ H + p reaction using microscopic nuclear reaction theory approach 



see 



19|)- 



Recently, it has been proposed to measure PV asymmetry of protons in n +^ He — t-^ 
H + p reaction with polarized neutrons at the Spallation Neutron Source at the Oak Ridge 
National Laboratory. For typical neutron energy En ~ 0.01 eV, which corresponds to a wave 
vector kn ~ 2.19 ■ 10~^/m~^, the energy of outgoing protons and proton wave vector are 
Ep = 0.764 Me and kp = 0.19 /m~^, correspondingly. Taking a characteristic ^He radius 
as R = 1.97 fm, one obtains {knR) ~ 4 ■ 10"^ and {kpR) ~ 0.4. Therefore, for the initial 
channel, contributions from p-wave neutrons to a reaction matrix (amplitude) are highly 
suppressed, whereas for the final channel, the amplitude with orbital momenta of protons 
I = and I = 1 has the same order of magnitude. The contribution from (i-wave protons 
is suppressed by a factor ~ 0.025; therefore, one can ignore d- waves within the accuracy of 
our estimates. 

PV asymmetry apy of outgoing protons, in directions along to neutron polarization and 
opposite to it, is proportional to PV correlation (a - kp). Using standard techniques (see, for 



example 20|), one can represent these asymmetries in terms of matrix R which is related to 



reaction matrix T and to S-matrix as 



R = 2mf =i-S. (8) 
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Then, for our case, 



apy = -Re\-3V2 < 01|i?^|10 > ■ < 00|i?°|00>* 



r 



+ (Ve < ll|i?°|00 > +6 < ll|i?^|10 >) < 10|i?^|10>*], (9) 



where 



r = (I < 00|i?°|00 > p + 3| < 10|i?^|10 > n. (10) 

We use spin-channel representation, where for the matrix element < s'l'\R'^\sl >, / and /' 
are orbital momenta of initial and final channels with corresponding spin-channels s and s' 
and J is the total spin of the system. Calculations of matrix elements < s'l'\R'^\sl > for PV 



effects in nuclear reactions have been done 



201 1 using distorted wave Born approximation in 



microscopic theory of nuclear reactions 2l|. They lead to the PV amplitudes induced by 
parity violating potential W, 

Rf;y = 2m<^^\W\^t>, (11) 



/ 



where "iffr are the eigenfunctions of the nuclear P-invariant Hamiltonian with the appropri- 



ate boundary conditions 



2l|: 



k 



Here, (f)k is the wave function of the k resonance and Xmi-^) is the potential scattering 
wave function in the channel m. The coefficient 

"'^-^^^''^ - (2vr)i E-E,±^, 

describes nuclear resonances contributions and the coefficient 6^^^ (i?, E') describes poten- 
tial scattering and interactions between the continuous spectrum and resonances. Here, E^, 
Ffc, and are the energy, the total width, and the partial width in the channel i of the fc-th 
resonance, E is the neutron energy, and 5i is the potential scattering phase in the channel z; 
(T l)^ = (27r)^ < Xi{E)\V\(j)k >, where V is a residual interaction operator. As was shown 



m [20| for nuclei with rather large atomic numbers, the resonance contribution is dominant. 
Then, for the simplest case with only two resonances with opposite parities, the expressions 
for matrix element R for neutron-proton reaction with parity violation is: 

^ ' ^" {E-E, + ^^,/2)(E - E, + ^r,/2) ' ' ^'^^ 



and with conservation of parity (one resonance contribution) is: 

< s'i'\R^\si >= ;(rr(^)rf.(sO)^ .(,.^,.) 

' ' {E-Ei + tTi/2) ^ ^ 

where w = — J (piW (piidr is PV nuclear matrix element mixing parities of two resonances. 

This technique has been proven to work very well for the calculation of nuclear PV effects 
for intermediate and heavy nuclei. Assuming the dominant resonance contribution to PV 
effects for the n +^ He -^^ H + p reaction, we use this approach to estimate characteris- 
tic values of PV effects using parametrization of PV effects in terms of known resonance 



structure of the system. Fortunately, the detailed structure of resonances {^He levels) [22 1 



231 ] are well known for this reaction from 



and low energy neutron scattering parameters 
numerous experiments. 

To estimate PV asymmetry in He — H + p reactions using the described formalism, 
we take into account all known resonances 



22 



231] which result in multi-resonance represen- 



tation for R matrix elements. From the selection rules for angular momenta (see Eq. (Q and 



general expressions in 



201]). one can see that for low energy neutrons only resonances with 



the total spin of J = 0, 1 contribute to PV asymmetries of the interest. Thus, we consider 



contributions from nine low energy resonances 22|, |23| (see Table ( IIIII) ): one resonance with 
total angular momentum and parity J'^ = 0"^, three with J'^ = 0~, four with J'" = 1~, 
and one with = 1+. For further calculations, we assume that all weak matrix elements, 
which mix resonances with opposite parities, have the same values and are described by a 



phenomenological formula 20] w = 2 ■ 10 "^eV ^/D{eV) (where D is an average energy level 
spacing). This formula is in good agreement with other statistical nuclear model estimates 



24j-l26| of nuclear weak matrix elements for medium and heavy nuclei. The extrapolation 
of this formula to the region of one-particle nuclear excitation leads to the correct value 
for weak nucleon-nucleon interaction. Therefore, one can use this approximation for rough 
estimates of average values of weak matrix elements in few-body systems. This leads to the 
value of weak matrix element w = 0.5 eV (with Z) ~ 6 MeV), which is rather close to the 
typical value of one particle weak matrix element. One can see from Eqs. fll4l) and (JT5l) that 
the expressions for PV and PC i? matrices depend not on neutron and proton partial widths 
but on their amplitudes, the values of which depend on particular spin-channels. Since we 
know only partial widths, we have to make assumptions about values of amplitudes of par- 
tial widths for a specific spin-channel and about their signs (phases). This leads to another 



uncertainty in our estimation in addition to the previously given assumption about weak 
matrix elements. To treat the spin-channel dependence of partial width amplitudes, we 
assume that partial widths for each spin-channel are equal to each other. This gives us an 
average factor of uncertainly of about 2. The signs of width amplitudes, as well as the signs 
of weak matrix elements w, are left undetermined (random). This also can lead to a factor 
of uncertainly of 2 or 3. Therefore, one can see that the uncertainly of our multi-resonance 
calculations is about of one order of magnitude. 



Taking into account these considerations and using resonance parameters 22|, |23| of the 
Table (IIIip . one can estimate the PV asymmetry for thermal neutrons as 



apy = -(1 - 4) ■ 10"^. (16) 
The set of resonance parameters of the Table OIVp results in a slightly lager PV asymmetry 

apv = -(4 - 8) ■ 10-^. (17) 



The difference between these two sets is related to the discrepancy between 23| and 22| for 
resonance parameters for the first positive resonance {En = 0.430 MeV). 

To show contributions of each resonance to PV asymmetry apy, we normalized contri- 
butions from each resonance in terms of relative intensity to the strongest one, which is 
taken as 100% (see last two columns in Tables ( IIII|) and (jIVP ). Some resonances contribute 
through two different spin-channels, s = and s = 1. In those cases, the contributions from 
two spin-channels can be either with the same sign or with the opposite sign, depending on 
unknown phases of amplitudes of partial widths and weak matrix elements (see, for exam- 
ple, resonance at 3.062 MeV in Table (IIIip ). As can be seen from these tables, different 
resonances contribute essentially differently to the value of PV violating effects. Moreover, 
different sets of resonance parameters can change weights of the resonances for a particu- 
lar asymmetry. For example, the lowest 0~-resonance contribution to the asymmetry apy 
appears to be 3% using the parameters of Table flllll) . while it would be the dominant one 
using the parameters of Table (lIVp . This is related to the fact that for the set of Table (lllip , 
the contribution of the 0~-resonance to the apy is suppressed by a factor of about 40 due to 
destructive interference between parity conserving and parity violating amplitudes. There- 
fore, the readability of this method can be essentially improved by increasing the accuracy 
in measurements of parameters of the most "important" resonances. 

10 
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FIG. 1. (Color online) Resonance enhancement of of the apv asymmetry (for the first set of 
parameters). 



It should be noted that the estimated value of the PV asymmetry apv at thermal energy 
(see Eqs. flTB]) and flTTl)) is surprisingly in very good agreement with exact calculations for 
zero energy neutrons 27| . This could be considered as an additional argument for reliability 
of the suggested resonance approach. Also, matching the estimated value of the observable 
parameter with exact calculations at low energy gives us the opportunity to predict PV 
effects in a wide range of neutron energies. 



III. CONCLUSIONS 

Considering a variety of methods of the calculation of PV effects in low energy few-body 
nuclear reactions, the EFT approach can be a solution for the discrepancy between DDH 
description of PV effects and experiments. Furthermore, the "true" EFT approach, which 
involves a solution of AGS-type of few body equations, can be the most promising method 
without introducing any additional model dependencies or parameters. 

At the same time, the use of all other approaches is still useful in some cases. In partic- 
ular, the simple nuclear reaction approach could be very useful for a preliminary study of 
new PV effects in few-body system, since it does not require heavy calculations and provides 
estimates for both PV and parity conserving effects for rather wide energy regions. 
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TABLE I. DDH PV coupling constants in units of 10 ^. Strong couplings are |^ = 13.9, |^ = 0.84, 
1^ = 20, Hp = 3.7, and k^j = 0. 4-paramter fit and 3-parameter fit uses the same /i^ and with 
DDH 'best'. 



DDH CouDhne 


; DDH 'best' 4-parameter fit [15] 
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TABLE II. Parameters and 

represented by qa by using 

%j = T? - Ti ■ Tj). 



operators of parity violating potentials. ttNN coupling g-j^NN can be 
Goldberger-Treiman relation, = gAmN/F-;^ with F-^ = 92.4 MeV. 
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TABLE III. Resonance parameters (Set 1). Here, Er is a resonance energy; T and J'^ are resonance 
isospin and the total resonance spin with parity; T and Tp are total and proton widths; T^, K 
and I are neutron width, reduced width, and angular momentum, correspondingly; and apy {% ) 
is normalized contribution of the resonance to apv- 



Er{MeV) 


J"" 


l 


T 


TJMeV) 

ft \ ^ ' / 




VJMeV) 


T{MeV) 


apv (% ) 


-0.211 


0+ 










954.4 


1.153 


1.153 




0.430 


0- 


1 





0.48 




0.05 


0.53 


3.1 


3.062 


1- 


1 


1 


2.76 




3.44 


6.20 


100±26 


3.672 


1- 


1 





2.87 




3.08 


6.10 


75±24 


4.702 


0- 


1 


1 


3.85 




4.12 


7.97 


20 


5.372 


1- 


1 


1 


6.14 




6.52 


12.66 


79±18 


7.732 


1+ 








4.66 




4.725 


9.89 




7.792 


1- 


1 





0.08 




0.07 


3.92 


2±1 


8.0G2 


0- 


1 





0.01 




0.01 


1.89 


11 



TABLE IV. Resonance parameters (Set 2). Here, Er is a resonance energy; T and J'^ are resonance 
isospin and the total resonance spin with parity; T and Tp are total and proton widths; r„, F'^ 
and I are neutron width, reduced width, and angular momentum, correspondingly; apy {% ) is 
normalized contribution of the resonance to apv- 



Er{MeV) 




l 


T 


TniMeV) 


Kiev) 


TpiMeV) 


r{MeV) 


apv i% ) 


-0.211 


0+ 










954.4 


1.153 


1.153 




0.430 


0- 


1 





0.20 




0.640 


0.84 


100 


3.062 


1- 


1 


1 


2.76 




3.44 


6.20 


82±27 


3.672 


1- 


1 





2.87 




3.08 


6.10 


62±20 


4.702 


0- 


1 


1 


3.85 




4.12 


7.97 


16 


5.372 


1- 


1 


1 


6.14 




6.52 


12.66 


65±15 


7.732 


1+ 








4.66 




4.725 


9.89 




7.792 


1- 


1 





0.08 




0.07 


3.92 


2±1 


8.062 


0- 


1 





0.01 




0.01 


4.89 


1 
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